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ON MATRICES WHOSE ELEMENTS ARE INTEGERS. 

By Oswald Veblen and Philip Franklin. 

Introduction. 

1. The purpose of this article is strictly expository. The aim is to 
set forth some of the theorems on matrices whose elements are integers. 
These theorems have applications in Analysis Situs* and the systematic 
treatment of them directly in terms of integers here given will no doubt 
be useful to students of that subject. While the closely allied algebraic 
theory is to be found in B6cher's Introduction to Higher Algebra, and the 
matter here given is to some extent discussed in Muth's Elementartheiler 
and in Scott and Mathews' Determinants, there is no readily accessible 
treatment of the subject from the point of view here adopted. 

2. The object of our study will be a matrix of a rows and /S columns: 



(1) ^ = 11 6/ 11 = 



6ll 61^ ... 61^ 
I ^ 8 

€2 62 • • • 62 



The elements of E are integers. The term " integer " here includes 
negative integers and zero; but we shall assume that at least one element 
is different from zero. 

Our definition of the product of two matrices || «,-' || and || t;,-'' || is: 

(2) li v^' II = II e/ II • li -n^' II , 
where 

(3) p/ = lfei*-W. 

The number of rows of the second matrix must be equal to the number of 
colunms of the first; and the product has as many rows as the first matrix 
and as many columns as the second. If the matrices are square, the 
product will be square, and the determinant of the product will be equal 
to the product of the determinants of the factors. 

The inverse of a square matrix A, of determinant imity, will be the 
matrix A~^ such that: 

A-^-A = A-A-i = 7, 

* Cf. O. Veblen, Cambridge Colloqmum Lectures on Analysis Situs. 
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2 OSWALD VEBLEN AND PHILIP FRANKLIN. 

where I denotes the identity matrix || 5/ |] , a square matrix with all the 
elements in the main diagonal + 1 and all the remaining elements zeros. 
The element o/ of A~^ will evidently be the cof actor of a/ in the deter- 
minant of A. A is restricted to be of determinant unity to insure the 
elements of the inverse matrix being integers. 

Elementary Transformations. 

3. Let us consider two types of transformations of E: 

(a) To replace each element of the rth row {e/) by the element {er' + qej) 
where q is either + 1 or — 1 and s =^ r. This operation is described as 
adding the sth row to the rth row or subtracting the sth row from the 
rth row. 

(fe) To add a column to or subtract it from another column. 

The operation (a) is equivalent to multiplying E on the left by a 
square matrix of a rows Ao = || o.-' |i in which all the elements are zeros 
except those of the main diagonal which are + 1, and Or' which is q. For 
the expressions given by (3) for the elements of the product, i.e., 

(4) Pi' = Xa*-€k' 

reduce to the single term e/ except when i = r; in which case thej' give 
the two terms: 

e/ + qej. 

That is, the operation (a) transforms E into Aq-E. 

In like manner, the operation (b) corresponds to multiplying E on the 
right by a square matrix of /3 rows Bo = 1! b/ || in which all the elements 
are zeros except those of the main diagonal which are 1 and b/ which is q. 

If the operation (a) be repeated n times, where » is a positive integer, 
the effect is an operation identical with (a) except that q is replaced by 
the integer w or — n. Correspondingly, the effect of multiplying the 
matrix Ao bj^ itself repeatedly is to change the element a,' to ± n. 

The inverse oi Ao ii a/ = ±1 is the same matrix except that the 
sign of a/ is changed. Hence the inverse of an operation of type (a) 
is an operation of the same type. The determinant of A o is + 1. 

Similar statements hold with regard to the operation (6) and the 
matrix Bo- 

4. The operation of interchanging two rows of a matrix and changing the 
signs of all the elements of one of them can be expressed as a sequence of 
operations of type (a). For if we add the rth row to the sth, then subtract 
the sth row of the resulting matrix from the rth, and finally add the rth 
row to the sth, the elements of the rth and sth rows (and the qth column) 
will be, successively: 

(fr«, e,«); (€r«, €r« + €.«); (- e.«, «r« + e.«); (-«.«,«,«); 
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and the resulting matrix will thus be that obtained by changing the signs 
of the elements of the sth row and then interchanging the rth and sth rows. 
In Uke manner, the operation of interchanging two columns and changing 
the signs of the elements of one of them is expressible as a sequence of opera- 
tions of type {b) . 

5. In place of our two fundamental operations (a) and (b) we might 
have restricted ourselves to the operations: 

(a') To add a row to, or subtract it from, an adjacent row. 

(6') To add a column to, or subtract it from, an adjacent column. 

We reduce the operation (a) to a sequence of operations (a') in the 
following manner. For definiteness, let us speak of the sth row as follow- 
ing the rth; in the reverse case " following " is to be replaced by " pre- 
ceding " in the argument. Add each row from the rth to the (s — l)th 
inclusive to the next following row, beginning with the (s — l)th. Then, 
in the resulting matrix subtract each row from the (r + l)th to the 
(s — l)th inclusive from the following row, beginning at the (r -|- l)th. 
Next add each row from the (r + l)th to the (s — 2)th inclusive to the 
following, beginning at the (s — 2)th. Finally subtract each row from 
the rth to the (s — 2)th from the following row, beginning at the rth. 

If the reader will write out the expressions for the elements of the 
matrix in a single column and the rows affected (from the rth to the sth), 
he will find that the resulting matrix only differs from our original matrix 
in having its rth row added to its sth (or subtracted from it) . To perform 
the inverse operation we need only repeat the process, subtracting the 
rth row from the (r + l)th in the second step, and adding it to the 
(r + l)th in the last step; the other operations remaining as before. 

As a similar argument holds for steps (6) and (b'), if we replace rows by 
columns, we conclude that the transformations built up from steps (a) and 
(b) are no more general than those built up from steps (a') and (&')• 

Determinant Factors. 

6. Consider the set of 7-rowed (0 < 7 ^ a, y ^ 0) determinants 
which can be formed from E by omitting a — 7 of the rows and — y 
of the columns in all possible ways. The highest common factor (H. C. 
F.) of such a set of determinants, if the determinants are not all zero, is 
denoted by Dy and is called the 7th determinant factor* of E. 

The determinant factors are unchanged when the matrix is operated on 
by transformations of type (a) or (6) . For, consider the effect of an opera- 
tion of type (a) which consists in adding the rth row to (or subtracting it 
from) the sth, on the 7-rowed determinants in question. All such de- 

* Cf. Scott and Mathews' Determinants, p. 76. 
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terminants which do not contain elements from the sth row are obviously 
unaffected, while those that contain elements from both the rth row and 
the sth are not affected because of an elementary theorem on determinants. 
The remaining 7-rowed determinants, as A„ which contain elements from 
the sth row and not from the rth, are converted into determinants of the 
form Ay zt Ay where Ay' is the 7-rowed determinant obtained from Ay 
by replacing the elements from the sth row by elements from the same 
columns of the matrix and from the rth row. 
The proof for operations of type (6) is sinailar. 

7. The foUomng theorem has an application in Analysis Situs: If a 
matrix E is such that each column either consists entirely of zeros or contains 
just two elements different from 0, one + 1 and the other — 1, all the de- 
terminant factors of the matrix are + 1 or — 1. 

The theorem follows immediately from the definition of a determinant 
factor, if we observe that anj' 7-rowed determinant formed hy striking 
out (a — 7) rows and (0 — 7) columns of the given matrix has either 
two, none or one element in each column different from zero. If no 
column is of the third type the detemainant is zero, since the sum of all 
the elements in each column is zero. If there is a column of the third 
type we evaluate the determinant with reference to such a column and 
then evaluate the minor with reference to a column with a single non-zero 
element in the minor, and so on. In this way we either arrive finally at 
± 1 for the value of the determinant, or else come to a minor with two 
or no non-zero elements in each column, in which case the determinant 
is zero. 

Reduction to Normal Form. 

8. Let us now consider a series of reductions of the matrix E which 
can be effected by transformations of types (a) and (b). If the first 
column consists entirely of zeros, add one of the other columns to it. Thus 
by a transformation of type (b) E is converted into a matrix Ei which has 
at least one non-zero element in the first colunon. If the first element 
of the first column is zero, add a row which contains a non-zero element 
in the first column to the first row. Thus by a transformation of type (a) 
El is converted into a matrix E^ for which the element of the first row and 
column is not zero. 

We shall now prove that if this non-zero element ci^ is not a factor 
of all the elements of the matrix, we can, by a series of transformations 
of types (a) and (b), replace it by a numerically smaller element different 
from zero. 

First, if one of the elements in the first column, «^r, is not divisible by 
€1^, upon adding the first row to (or subtracting it from) the rth a number of 
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times equal to the largest integer in the quotient of «^ by €iS an element 
numerically smaller than ei^ is obtained in the first column and rth row. 
Then, on subtracting the rth row from (or adding it to) the first row the 
matrix is converted into one with a smaller non-zero element in place of 
ei^. This has been done by a succession of operations of type (a). Sim- 
ilarly, if there were an element in the first row which did not contain ei^ 
as a factor, transformations of type (b), strictly analogous to those of 
type (a) just described, could be set up which would reduce the numerical 
value of ei^ 

Second, if «/ is a factor of all the elements of the first row and first 
column, but is not a factor of the element in the rth row and sth column, 
tr", we proceed as follows. Upon subtracting the first column from (or 
adding it to) the sth ei'/«i^ times (transformations of type (6)), the first 
element in the sth column becomes zero, while the rth is still not divisible 
by €iS since it has been changed by a multiple of ei^. If we now add the 
sth column to the first (an operation of type (6)), the element in the first 
row and column remains €i\ while the rth element in the first column is 
now not divisible by ei^ Hence we may replace €i^ by a numerically 
smaller element by the method of the preceding paragraph. 

If the element which replaces ei^ is not a factor of all the elements of 
the matrix, it may be still further reduced by a repetition of the process 
described in the two paragraphs above. If this process be continued, 
we must arrive after a finite number of steps — the number being less 
than the absolute value of ei^ — at a matrix whose first element di is a 
factor of all the others. When this point is reached, we may reduce all 
the elements in the first column except the first to zeros by operations of 
type (a), for we have merely to add the first row to (or subtract it from) 
any other row the number of times the first element of this row contains di. 
The elements of the first row, with the exception of the first, may be re- 
duced to zeros by similar operations of type (b). It is evident that all the 
elements of the matrix thus obtained contain the first elem ent as a factor. 

Thus we arrive at a matrix Ei in which the first element di of the 
first column is the H. C. F. of all the elements of E3 and in which all the 
other elements of the first row and of the first column are zero. By § 6, 
di is the H._C. F. of all the elements of E. 

9. Let Ez be the matrix obtained from E3 by deleting its first row and 
fijst column. By § 8, E3 may be reduced to a matrix with a leading ele- 
ment which is the H. C. F. of all its elements, and having all the other 
elements of the first row and column zero. 

As the transformations of types (a) and (b) which effect this reduction 
on E3 determine transformations of E3 of the same type, which leave its 
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first row and first column unchanged, we maj' reduce the matrix Ez to 
a matrix E4 in which the first element of the main diagonal, di, is the 
H. C. r. of all the elements of the matrix, the second element of the main 
diagonal, ^2, is the H. C. F. of all the elements except di, and all the remain- 
ing elements of the first two rows and first two columns are zero. 

B3' a continuation of this process we arrive bj^ a finite sequence of 
operations of types (o) and (6) at a matrix: 



(5) 



E* = 



\di 
di 
















dr 
















in which all the elements are zero except a sequence of elements di, 
(0 < i :^ r) common to the ith row and column and such that di is the 
H. C. F. of all d/s such that i^j^r. 

The di's may be positive or negative integers. We can make all ex- 
cept the last positive by a sequence of operations of type (a). For if 
di is negative, and we interchange the ith and rth rows, changing the sign 
of the elements in the rth, a permissible operation by § 4, and repeat the 
process, we arrive at a form in which di is positive and dr has changed 
sign. We may thus obtain a form in which d, {i < r) is positive, and dr 
will be positive or negative according as the number of negative signs in 
the form we started with was even or odd. We shall take this matrix, 
with at most one negative element, as the normal form E* in the dis- 
cussion which follows. 

Each operation of tj^pe (a) amounts, according to § 3, to multipljdng 
the matrix to which it is appUed on the left by a square matrix of type Aq 
of a rows, and each operation of type (6) amounts to multipljdng the 
matrix to which it is applied on the right by a square matrix of type Bo 
of 13 rows. Hence 
(6) E* = AE-B, 

where J. is a product of matrices of type Ao and B a product of matrices 
of type Bq. It is to be noted that the determinants of A and B are each 

+ 1. 

Let us introduce the notation Z>i = di, D2 = di-d^, ■ ■ • Dr = di-di- 
• ■ ■ dr, and observe that D^ (0 < 7 :^ r) is the H. C. F. of all the 7-rowed 
determinants which can be formed by striking out a — y rows and 
/3 — 7 columns from jE^*. That is, referring to § 6, they are the successive 
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determinant factors of E*. Since E* was derived from E by operations 
of types (o) and (6), they are also the determinant factors of E. 

Since the D/s are invariant under transformations of types (a) and 
(b), the di's, which are the quotients of successive D/s, (di+i = Di+i/Di), 
are also invariant under these transformations. They are called in- 
variant factors or elementary divisors.* 

The number r is also invariant under transformations of types (a) and 
(b) and is called the rank of the matrix E. 

The Matrices of Transformation. 

10. In the special case where £ is a square matrix of a rows whose 
determinant is + 1, equation (6) implies that the determinant of E* is 
+ 1. Hence r = a, and the numbers d.- must be + 1. 

We therefore have: 

(7) A-EB=I, 

in which A is a product of matrices of type Ao, B a product of matrices 
of type Bo and I is the identity matrix. We may write (7) in the form 

(8) E = A-^-I-B-^ = A-'-B-K 

It is evident from § 2 that when a = (3 every matrix of type Bo can be 
regarded also as one of type Ao, and the same is true of matrices inverse 
to those of types Ao or Bo. As the above equation shows that E is equal 
to a product of such matrices, we have the theorem: Any square matrix 
of determinant unity is expressible as a product of matrices which may be 
considered to be of type Ao, or to be of type Bo. 

Hence to multiply a matrix E of a rows and /3 columns on the left by a 
square matrix of a rows and determinant unity is equivalent to operating on 
E by a sequence of operations of type (a) ; and to multiply E on the right by 
a square matrix of /3 columns and determinant unity is equivalent to operating 
on E by a sequence of operations of type (b) . 

Also, since we may write (8) in either of the forms: 

(9) B-A-E = I or E-B-A = /, 

it follows that if E is a square matrix of determinant unity, it may be reduced 
to the form I by operations on rows only, or by operations on columns only. 

11. In the case of a general matrix E, we have from (6) 

(10) A-E = E*-B-\ 

Since the determinant of B~^ is 1, the H. C. F. of the elements of its first 

* We shall use the term invariant factor, following Bdcher, Introduction to Higher Algebra, 
pp. 269-70, since the term elementary divisor is sometimes used in another sense. 
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row is 1. Hence the H. C. F. of the elements of the first row of the 
matrix E*-B~^ is d\. As a similar statement applies to the remaining 
rows, we have the theorem: 

The matrix A has the property that the H. C. F. of the elements of the rth 
row of the matrix A-E is dr, the rth invariant factor of E. 

This suggests a method of building up A by means of the theorems: 

(1) That for any set of integers ey^ {0 < j :S a), a set of integers 
ffli^ (0 < j — a) can be found which are relatively prime and such that 

J^ai'-e/ = ti 

where ti is the H. C. F. of the a e/'s; and 

(2) That there exists a matrix A of determinant unity with the numbers 
Ui' as the elements of its first row. 

The derivation of equation (6) by this method is longer than that 
given in §§ 8 and 9 and is therefore omitted. 

Diophantine Equations of the First Degree. 

12. Consider the problem of finding the integral solutions of the 
f ollomng set of equations : 

ti^Xi + ei^a^o + • • • + «i% = Pi, 

QjN (i^Xi + ii^Xn + • • • + «2% = p2, 

^a^l + ea^ajo + • • • + t^^Xp = p„. 

If X denotes the matrix of one column, and jS rows 






x$ 



and P a similar matrix with pi, Pi • • • p^ as the elements of its one column 
and a rows, equations (11) may be written 



{i-Z) 


ii-A = f. 


But from (6) we have 




(13) 


E = A-^-E*-B-\ 


and consequently 




(14) 


A-^-E*-B-^-X = P, 


or 




(15) 


E*-B-^-T= AP. 
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Let us set Q = A-P, a, matrix of one column and a rows, and denote 
its elements by qi, Qi ■ • ■ q^. Also let yi, 2/2 • • • t/^j be the elements of the 
matrix Y = B~^-X, which is of one column and ^ rows. Then (15) 
becomes : 
(16) E*-Y = Q, 

which is equivalent to the set of a equations: 

(17) 



diVi = qi 
= gy 



(0 < i ^ r), 
(r < i < a), 



where r is the rank of E. If equations (17) are to be consistent, the q/s 
must all be zero, and in this case the solution is. 



(18) 



Vi = 



_ Qi 



di 



yj is arbitrary 



(0 < i < r), 
{r <j^ /3). 



To express the condition that equations (17) be consistent and solvable 
in integers, in terms of the coefficients of (11), we proceed as follows. 
Form the " augmented matrix " of the system, a matrix E oi a rows 
and |3 + 1 columns whose ith row has as its elements : 



«< , «<■ 



•, e." 



Vi- 



The matrix S formed by multiplying E on the left by A will have as the 
elements of its iih. row {Q < i -^ a): 



, s, , qi, 



Si , Si , 

where the s/'s are the elements of the matrix: 

S = lis/- II = A-E. 

Since multiplying the matrix S by B reduces it to the normal form, it may 
be reduced to the form E* by operations on columns only; which shows 
that S may be reduced by operations on columns only to the form: 



(19) 



E* = 



d^ 





di 

















dr 



-gi 
-32 



-qr\ 



••• ••• -g., 
In order that (11) be solvable at all, we found that qi must be zero for 
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values of i greater than r. This shows that the rank of £■* is r. If in 
addition we require the solutions to be integers, g, must be divisible by 
di for i < r. Hence E* may be reduced to normal form bj^ adding the 
ith column to the last g,/di times. Hence its invariant factors must be 
the same as the d/s, i.e., those of E. Conversely, if this condition is 
satisfied, each q, will be divisible by the corresponding d,-, and the solu- 
tions of (11) will be integers. _ 

Since E* was obtained from E by elenaentarj^ transformations, it 
has the same rank and invariant factors as E. Hence we have proved 
the two theorems: 

A necessary and sufficient condition that the equations (11) have a set of 
integral solutions is that the augmented matrix E have the same rank and 
invariant factors as the matrix of the coefficients E. 

13. Since the solutions of (17) are given by (18), and since X = BY, 
the solutions of (11) are: 

(20) Xi = if b/?/y = £ 5/ U + Jf 5/t/y (0 < i ^ /?), 

j=i j=\ dj j=>-+i 

in which yr+u Vr+i, • • ■, ye S'^e arbitrary integers. 

If the equations were homogeneous, the p,'s would all be zero, and 
hence the g/s would also be zero. Hence the solutions would be of the 
form: 

(21) Xi = jfbJyy (0 < i ^ 0), 

in which yr+i, yr+2, • • ■, yg are arbitrary integers. 

Consequently, for such equations we have the theorem: 
A set of linear homogeneous equations whose coefficients are integers has 
a set of j3 — r linearly independent solutions each of which is a set of relatively 
prime integers, if is the number of unknowns and r the rank of the matrix 
of the coefficients. All other solutions in integers are linearly dependent on 
these — r linearly independent solutions, the coefficients of the linear rela- 
tions being integers. 

This result was to be expected, since if a set of linear homogeneous 

equations are solvable in rational numbers, they are solvable in integers. 

By comparing (20) and (21) we obtain the further result: 

// one set of integers satisfying equations (11) be given, the other solutions 

are obtained by adding to it the solutions of the homogeneous equations which 

result when the right members of (11) are replaced by zeros. 

The theorems of this paragraph were first given in complete form by 
H. J. S. Smith,* although he was anticipated to some extent by Heger.f 

* Smith, H. J. S., On Systems of Linear Indeterminate Equations and Congruences, Philos. 
Transactions, Vol. 151, pp. 293 f. Collected Works, XII, pp. 367 ff. 

t Heger, Ignaz, Mem. Vienna Academy, Vol. XIV, second part, p. 111. 



ON MATRICES WHOSE ELEMENTS ARE INTEGERS. 



11 



Skew-Sjrattmetric Matrices. 
14. A skew-symmetric matrix is one in which 

(22) 

Let us define as the conjugate of a square matrix the matrix obtained from 
it by interchanging rows and columns. Evidently if a skew-symmetric 
matrix be pre-multiplied by any square matrix, and post-multiplied by 
the conjugate of this matrix, it will remain skew-symmetric. 
If A is the matrix defined in § 9 such that 

(6) A-E-B = E*, 

the matrix A-E, by § 11, has di as the H. C. F. of the elements of the 
first row. Since multiplication on the right corresponds to operations on 
columns only and leaves the H. C. F. of the elements of the first row un- 
changed, di is also the H. C. F. of the elements of the first row of A -E-A' 
where A' denotes the conjugate of A. Since A-E- A' is skew-symmetric, 
di will also be the H. C. F. of the elements of its first column. 

We reduce A-E -A' further as follows: If the second element of the 
first row does not divide all the remaining elements in that row, let ei^ be 
one which it does not divide. Subtract the second column from (or add 
it to) the ith a number of times equal to the greatest integer in the quotient 
ei-'/er, thus replacing «f' by an element numerically less than er. Upon 
subtracting the jth column from (or adding it to) the second, we obtain 
an element in the first place of the second column smaller than the one 
there before. xlU these operations leave the first column unchanged, and 
since the matrix was skew-symmetric, a similar set of operations on the 
rows reduces the matrix to a skew-symmetric matrix with the first element 
in the second column numerically smaller than before. By repeating 
these operations a sufficient number of times — at most | ei^ | times — this 
first element will be the highest common factor of the elements in the 
first row, and consequently of the elements of the matrix. When this 
condition is reached, we combine the second column with the other 
columns such a number of times that all the elements in the first row after 
the second will be zero, and perform similar operations on the rows. Then 
we combine the first column with the other columns such a number of 
times that all the elements in the second row after the first wiU be zero, 
and perform similar operations on the rows. This will reduce our matrix 
to the form: 



(23) 



E^ 






dx 








-d. 




















ts' 








€/ 













«4'' 
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The matrix obtained from Ei by deleting its first two rows and columns 
is skew-symmetric, and by apphang the above process to it, in a way 
wholly analogous to the way bj^ which we extended our initial process of 
reduction in § 9, we may, by a finite number of operations, reduce our 
matrix to one, E = \\ €»-"' |], in which 

ezi-i^' = di-, €2i2'-i = - di (0 < i ^ p) 

and the remaining elements are zero. That is, our matrix consists of a 
series of skew blocks of two non-zero elements each along the main 
diagonal, surrounded by zero elements. 

Since in the above process we have always performed identical opera- 
tions on rows and columns, we may write : 

(24) E = V-E-U' 

where U and U' are conjugate matrices whose determinants are -f 1. 

Since interchanging the first and second, third and fourth, • • • (2n 
— l)th and 2nth rows, and changing the signs of the even rows would 
reduce this matrix to the usual normal form E*, the d/s appearing in E 
must be identical with those of E*, i.e., the invariant factors of E. Hence, 
we have the result: 

The invariant factors of a skew-symmetric matrix are equal in pairs, and 
the rank of such a matrix is an even number. A skew-symmetric matrix 
may be reduced to the " skew " normal form, E, by multiplying on the left 
by a unimodular matrix U and on the right by its conjugate, U'. 

Symmetric Matrices. 
15. A symmetric matrix is one in which: 

(25) «/ = e/. 

Since a symmetric matrix retains its symmetry when we perform any 
operations on its rows, provided we perform the same operations on its 
columns, the question naturally arises whether a process similar to that 
of the preceding paragraph exists which will enable us to reduce such 
matrices to their normal form by means of a matrix and its conjugate. 
This question must be answered in the negative,* as is proved by the 
following example. The matrix 

i;2 1! 

:i 21 



(26) E = 



' H ' 



* On p. 189 of Scott and Mathews' Determinants the erroneous statement is made that 
symmetric matrices with integral elements can alwaj's be reduced to normal form by identical 
operations on rows and columns. 



E* = 


1 




3 


u = \ 


a 
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can not be reduced to its normal form, 

(27) 

by a matrix 

(28) 

and its conjugate U', where a, b, c and d are integers, since one of the 
conditions the elements of U would have to satisfy is: 

(29) 2a- + 2ab + 26^ = 1. 

The reduction and classification of symmetric matrices by identical 
operations on rows and columns is thus a problem of a different order 
from those which have been considered in this paper. Equivalence under 
such operations involves much more than the equivalence of invariant 
factors. This classification is the fundamental problem of the arithmetic 
theory of quadratic forms.* 

Matrices with Elements Reduced, Modulo 2. 

16. In many applications of matrices to Analysis Situs, it is found 
convenient to reduce the elements of the matrices modulo 2. On re- 
ducing modulo 2, the equation (6) becomes 

(30) A-EB = E*, 

in which E can represent an arbitrary matrix of a rows and /S columns 
whose elements are and 1, A and B represent square matrices of de- 
terminant unity (mod. 2) of a and /3 row^s respectively, and E is a matrix 
all of whose elements are except a sequence of elements along the main 
diagonal which are 1. This follows from the fact that if one of the in- 
variant factors of E* is even, so are all the following invariant factors 
since they contain this one as a factor. The number of I's in E* is the 
rank of E*. It is less than or equal to the rank of E*, and differs from 
it by the number of even invariant factors of E. 

Ssmunetric Matrices, Modulo 2. 

17. The theory of symmetric matrices, mod. 2, is not subject to the 
difficidties referred to in § 15. The reduction of such a matrix to normal 
form may be effected as follows: First interchange rows (performing the 
same interchange of columns) until the main diagonal consists of a series 
of I's followed by a series of O's. This can be effected by elementary 

* Cf. Encyclop6die des Sciences Math^matiques, Tome I, vol. 3, p. 101. 
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transformations according to § 4 because the negative of any element is 
the same as the element itself, modulo 2. Add the first row to every row 
whose first element is a 1, and the first column to the corresponding 
columns. Repeat this for the second row and column performing a new 
interchange of rows and columns, if necessary, and continue until there 
are no elements different from zero in the main diagonal after those used. 
The part of the matrix still to be normalized is now in the skew- sym- 
metric form, since + 1 = — 1 (mod. 2) and may be normalized by the 
process of § 14. Thus bj^ identical operations on rows and columns we 
have reduced our matrix to the form E* = || ej \\ in which: 

€.■•• = 1 {0 <i^p); ep+2,-i''+2' = 1; €^21"+^'-' = 1 (0 < i < <?), 

and all the remaining elements of the matrix are zero. That is, the non- 
zero elements consist of a series of I's in the main diagonal, followed bj^ 
a series of skew blocks, each containing two I's. If p = 0, this noatrix 
can not be reduced further; but if p + 0, it maj' be reduced to a form con- 
taining one or two I's in the main diagonal (according as p is odd or even) 
and a series of skew blocks, or to a form containing a series of I's in the 
main diagonal and no skew blocks. This further reduction depends on 
the fact that a group of three I's in the main diagonal of a matrix in the 
above form may be replaced by a single 1 in the main diagonal and a skew 
block of two. 

The steps of the process in the case of a three-rowed square matrix 
are, first adding the first row and column to the second row and column 
respectively, then adding the third row and column to the first row and 
column respectivel}', and finally adding the second row and column to 
the third. The matrix becomes successivelj^ : 



1 





1 








1 



1 


1 





1 














1 



! 1 1 
I 1 
I 1 1 



io 1 0: 
!i 

io li 



These steps may be made in reverse order to effect the inverse transforma- 
tion, and are obviously typical of the steps which can be applied to any 
matrix E* for which p > 0- 

Symmetric Matrices, Modulo p. 

18. In the case of an odd prime modulus, the reduction of a sym- 
metric matrix is even simpler than in the modulo 2 case, and the normal 
form is a matrix all of whose elements are except a sequence down the 
main diagonal. For, by a set of interchanges on rows, followed by 
similar ones on columns, we can obtain a non-zero element in the first row 
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of the matrix. If this is not the leading element, by adding the column 
containing it to the first column, and the corresponding row to the first 
row, we obtain a leading element, ei\ which is not zero. We may then 
add the first row to the others so as to make all the other elements in the 
first column 0, and operate similarly on the columns. The number of 
times, n, we must add the first row to a row with first element tr^ is given 
by solving the congruence: 

(31) neii + e.i = (mod. p) 

which has a root since p is prime. The reduction is continued as in § 9. 
By interchanges of rows and columns after we have reduced our matrix 
to a form similar to that given by (5), we can change the order of the d/s, 
and have to decide on some definite order to get a single normal form. 
We may, for example, write the p — 1 non-zero elements of our system as 
the integers from 1 to (p — l)/2 with plus and minus signs; and take as 
the normal order that of absolute value of the elements when written in 
this form. 



